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ABSTRACT 


V 


The structure of premixed flames propagating in combustible systems, 
containing uniformly distributed volatile fuel particles, in an oxidizing gas 
mixture is analyzed. This analysis is motivated by experiments conducted at 
NASA Lewis Research Center on the structure of flames propagating in 
combustible mixtures of lycopodium particles and air. Several interesting 
modes of flame propagation were observed in these experiments depending on 
the number density and the initial size of the fuel particle. The experimental 
results show that steady flame propagation occurs even if the initial 
equivalence ratio of the combustible mixture based on the gaseous fuel 
available in the particles, <|) u is substantially larger than unity. A model is 

developed to explain these experimental observations. In the model it is 
presumed that the fuel particles vaporize first to yield a gaseous fuel of known 
chemical composition which then reacts with oxygen in a one-step overall 
process. The activation energy of the chemical reaction is presumed to be 

large. The activation energy characterizing the kinetics of vaporization is 
also presumed to be large. The equations governing the structure of the flame 
were integrated numerically. It is shown that the interplay of vaporization 
kinetics and oxidation process, can result in steady flame propagation in 
combustible mixtures where the value of <J) U is substantially larger than unity. 

This prediction is in agreement with experimental observations. 

L Introduction. 

The aim of this study is to provide a numerical and asymptotic 
description of the structure of planar laminar flames, propagating in t 
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medium containing a uniform cloud of fuel-particles premixed with air. 
Attention is restricted here to systems where the fuel-particles first vaporize 
to form a known gaseous fuel, which is then oxidized in the gas-phase. This 
research was supported for the period September 14, 1991 to September 13, 
1992. This research was performed in collaboration with Professor A. Linan at 
Instituto Nacional de Tecnica Aeroespacial and Escuela Tecnica Superior de 
Ingenieros Aeronauticos, Madrid Spain. 

The model developed here is built on previous asymptotic analysis 
structure of premixed particle-cloud flames [1]. In this previous analysis the 
vaporization rate of the fuel-particles was presumed to be proportional to 

T s 133 , where T s is the temperature of the fuel-particles. Hence, the activation 
temperature in an Arrhenius-type rate law for the vaporization process is not 
large. For simplicity, the temperature of the fuel-particles was presumed to be 
equal to the gas temperature. The analysis was performed for values of <t> u > 1.0, 
where <|) u is the equivalence ratio based on the gaseous fuel available in the 
particles. For given values of (|) u the analysis yields results for the burning 

velocity and <J) g , where <(>g is the effective equivalence ratio in the reaction 
zone [1], The analysis shows that even though <J) U > 1.0, for certain cases the 
calculated value of <()g is less than unity. 

A two-temperature model is described here, where the temperature of 
the fuel-particles is different from the temperature of the gas. However the 
size of the fuel-particle is assumed to be small enough so that its velocity is 
nearly equal to the gas velocity. A general treatment of flame propagation 
supported by volatile fuel-particles should consider both radiative and 

molecular transport mechanisms, the temperature difference and the 

differences in velocity between the gas and the particulates, the 

nonadiabaticity of combustible systems of finite size, as well as the detailed 

kinetics of oxidation and the kinetics of vaporization pyrolysis. This work has 
the somewhat limited aim of examining the interplay of vaporization kinetic 
and oxidative kinetic processes for cases where nonadiabatic and radiative 
transport mechanisms are not substantial. 
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2. Formulation. 

A model is developed to describe steady, one-dimensional, planar flame 
propagation in a combustible mixture consisting of uniformly distributed fuel- 
particles in air. The initial number density of the particles, n u (number of 
particles per unit volume) and the initial radius a u are presumed to be known. 
All external forces including gravitational forces effects are assumed to be 
negligible. The kinetics of vaporization of a fuel-particle is represented by 
the expression 

m = 47ia 2 Aexp|-^j = - (!) 

where rh represents the mass of gaseous fuel vaporized per unit time from the 
fuel-particle, a and p s are the instantaneous radius and density respectively of 
a fuel-particle. The rate parameters A and T av which respectively represent 
the frequency factor and the activation temperature of the vaporization 
process are presumed to be known. The ratio T av /T s is presumed to be a large 
quantity. 

Consider a vaporizing fuel-particle in an ambient atmosphere. The 
balance equation for the mass fraction of the gaseous fuel vaporizing from the 
fuel-particle Yp, and the energy balance equation can be written as, 

__liL_^XE.X j_( r 2 k ^Xf\ = (2) 

4 k r 2 dr r 2 dr \ C P / 

m (3) 

4 n r 2 dr r 2 dr ' ° r ' 

where r is the radial coordinate. For the gas mixture k is the thermal 
conductivity, c p the heat capacity, and T the temperature. The quantity Lf = 
k/(pc p Dp) is the Lewis number of the gaseous fuel which is presumed to be a 
constant and Dp is the diffusion coefficient of the gaseous fuel. Boundary 
conditions for Eqs. 2 and 3 far from the surface of the fuel-particle, and at the 
surface of the fuel-particle can be written as. 
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T = T„, Y f = 
T = T S , m(l - Yp) = - 


Y F oo as r — 

4 n k t 

Cp L F 9r 


r = a. 


(4) 


The values of an< ^ ^F,~ will be determined from the solution of the 
equations governing the structure of the laminar flame. The temperature of 
the fuel-particle T s can be calculated from the energy balance equation 

q g = |4 k r 2 k = rh Ly + j it a 3 p s c s (5) 


where t represents the time, c s represents the heat capacity of fuel-particle, 
and L v represents the heat required to vaporize unit mass of vapor from a fuel- 
particle. The values of c p and c s are presumed to be constant. If the value of k 
is also presumed to be constant, then from integration of Eqs. 2 and 3, using Eq. 
4 to determine the constants of integration, the profiles of Yp and T can be 
expressed as 


and 


Y F = 1 -(1 - Y F co )exp 


m c p L f 
4 7C k r I’ 


T = Tc - 


c p rh 


+ IT„-T S + 


q g 

c p rh 


exp 


mc p 
4 7t k r 


Introducing the definition, 


x= (r L ^ _ ) > 

14 7C k r If = a 


( 6 ) 

(7) 


( 8 ) 


from Eqs. 5 and 7 it follows that 


4 tc a k X. (T^ -Ts) _ 47 iakLyA. 
exp (x) - 1 Cp 



it a 3 Ps Cs 


dTs 
dt ’ 


(9) 
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The term on the left side of Eq. 9 denotes the rate of heat transfer from 
the gas to the fuel-particle, with the effective Nusselt's number denoted by 
A, /[exp (A.) - 1], The first term following the equality sign on the right side of 

Eq. 9 denotes the latent heat required to vaporize the fuel-particle, and the 
second term denotes the sensible heat required to raise the temperature of the 
fuel-particle. It follows from Eqs. 1, and 8 that 

^Aexp(-I^). (10) 

Let T v denote a characteristic value of the temperature of the fuel- 
particle for which the value of A evaluated from Eq. 10 is of order unity. For T s 
< T v , the vaporization rate and consequently the value of A, is exponentially 
small. In the limit of small values of A, A/[exp(A) - 1] = 1. Hence the first term 
on the right side of Eq. 9 can be neglected in calculating the value of T s . For T s 
= T v , A is of order unity. If the value of T s is allowed to increase beyond T v , 
then the value of A would become exponentially large and the term on the left 
side of Eq. 9 would become exponentially small, consequently the rate of heat 
transfer from the gas-phase to the solid, will be small. Therefore, for values of 
T s beyond T v , the value of the derivative dT s /dt will be negligibly small, so that 
the value of A remains of order unity. For simplicity, it is presumed that rh = 0 
for T s < T v and the value of T s does not increase after it attains the value T v . 

The governing equations for the structure of the flame can be written 
as. 


V = p u u u + 1 7t p s n u u u (a 3 - a 3 ), 

d(V Yp) k u d 2 Y F _ n uPu u u ^ p„ 
dx Cp Lp dx 2 V p 

d(VYa) wn pu 

dx CpLo, dx* ’ P’ 


^^. ko d2X = ^ qp+QRWF ^ 

dx dx 2 V p 


(ID 

( 12 ) 

(13) 


( 14 ) 



Here the subscript u denotes conditions in the ambient reactant stream. 
The subscript 00 has been dropped from the quantities Yp, Y02 and T appearing 
in Eqs. 12-14. The quantities wp and wo 2 denote the rate of consumption of 
gaseous fuel and oxygen per unit volume, and Qr is the heat released per unit 
mass of gaseous fuel consumed. The independent variable x is related to the 

spatial coordinate x' as, 


x = 



dx'. 


(15) 


In the regions where T s < T v , the vaporization rates are assumed to be 
negligibly small hence the radius of the fuel-particle does not change and its 
temperature T s and the quantity q p can be calculated from the relations, 


and 


Ps c s V dT s 
Pu 


3 k U (r p2 -p2\ 

2T U ' s/ ’ 



Ps C s 


V dT s 

Pu dx 


(16) 

(17) 


When the value of the temperature of the fuel-particle reaches T v , they 
begin to vaporize and the increase in their temperature is negligibly small. 
The radius of the fuel-particle and the quantity q p can be calculated from the 

expressions. 


and. 


a da. = _ 

dx 


ku Pu 

Tu Cp Ps V . 


T - T v + 




Cjj - Ty) ) 
Lv /J ’ 


q p = rh (c p T v - L v ), m = 


P 



V d(a 3 ) 
p u dx 


(18) 

(19) 


In deducing the system of equations 12-19 the product pk was presumed 
to be a constant. The chemical reaction between the gaseous fuel and the 
oxidizer is presumed to occur by a one-step process. The system of equations 
12-19 were solved numerically. 


3. Results 


For purpose of illustration calculations were performed assuming that 
the gaseous fuel that evolves from the fuel-particle is methane. The chemical 
kinetic rate parameters characterizing the gas-phase oxidation were chosen 
such that the burning velocity for stoichiometric mixture of gaseous methane 
and air is 37 cm/s. The values of some of the parameters used in the 

calculations are k u = 0.00005 cal/(cm s K), p u = 0.001135 g/cm 3 , p s = 1.0 g/cm 3 - 
T u = 298 K, T v = 600 K, L v /Qr = 0.01. Results for the burning velocity u u as a 
function of the equivalence ratio ([> u , calculated based on the gaseous fuel 

available in the fuel-particles, for values of a u equal to 1.0 pm, 3.0 pm, and 5.0 

pm are shown in Fig. 1. Figure 1 shows that steady flame propagation is 

possible in fuel-rich combustible mixtures and the value of u u initially 

increases with increasing values of <]) u attains a maximum value then begins to 
decrease. In Fig. 2 the structure of the flame is plotted at a value of <|) u = 1.8 and 
a u = 3.0 pm and corresponds to the region where the value of u u increases with 
increasing values of <j> u - The various non-dimensional quantities appearing in 
this figure are defined as, 


_ Pu U u Cp _ T - T u q _ T s - T u 
k u 1 b ' * u 

_JL 


a = 


yf = ’ y°^ = ^ 

ICbU *02U 


( 20 ) 


where Tb is the adiabatic flame temperature. Figure 2 shows that gaseous fuel 

is the limiting reactant in this region. In Fig. 3 the structure of the flame is 
plotted at a value of <]> u = 3.69 and a u = 3.0 pm and corresponds to the region 
where the value of u u decreases with increasing values of <t> u . Figure 3 shows 
that oxygen is the limiting reactant in this region. 


4. Future Research. 


Future research will be directed toward performing numerical 
calculations using the system of equations 12-18 over a wide parametric range. 
Asymptotic description of the flame structure will also be attempted. 
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Figure 1. The burning velocity u u (cm/s) as a function of the 
equivalence ratio <{>u, for various values for the initial radius of the 
fuel-particle a u . 
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Figure 2. The structure of the flame for <]> u = 1.8 and a u = 3.0 






Figure 3. The structure of the flame for <J>u = 3.69 and a u = 3.0 |im. 



